Almost-stationary motions and gauge conditions in General Relativity 
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An almost-stationary gauge condition is proposed with a view to Numerical Relativity applica- 
tions. The time lines are defined as the integral curves of the timelike solutions of the harmonic 
almost-Killing equation. This vector equation is derived from a variational principle, by minimizing 
the deviations from isometry. The corresponding almost-stationary gauge condition allows us to put 
the field equations in hyperbolic form, both in the free-evolution ADM and in the Z4 formalisms. 
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I. INTRODUCTION 

It is well known that there is no preferred set of coor- 
dinate frames in Einstein's theory of gravitation, which 
is known as General Relativity precisely because of this 
fact. In practical applications, however, one is forced to 
consider specific values of the gravitational field compo- 
nents (spacetime metric, curvature tensor) and this can 
be done only after choosing a specific coordinate system. 
One just expects that this gauge choice will not hide the 
physics of the problem behind a mask of non-trivial co- 
ordinate effects. 

Harmonic coordinate systems have deserved much in- 
terest since the very beginning of Einstein's theory 0, ■ 
The coordinates themselves are defined by a set of four 
harmonic spacetime functions { $ M } , that is 



□ ^ = 



(1) 



where the box stands for the wave operator acting on 
functions. In an harmonic coordinate system one takes 



so that Einstein's field equations 



8 7T T„„ 
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get a very convenient form, in which the principal part 
is simplya scalar wave equation for every metric compo- 
nent 0, Q • This fact has been used by Yvonne Choquet- 
Bruhat for proving the well-posedness of the Cauchy 
problem in General Relativity [3J. 

In the recent years, this simplification is being used 
for building up Numerical Relativity codes, with in- 
teresting related developments also on the theoretical 
side 0,1a, IS O n the applications side, the 'gauge 
sources' variant of the harmonic gauge @ has been used 
by Pretorius in a numerical simulation of the evolution 
of a binary black-hole system: a full quasi-circular orbit 
has been achieved [UEy]- In this simulation, however, 
a grid velocity has been introduced in order to get into 
the co-rotating frame. This means that computational 
nodes are rotating with respect to the coordinate system. 
In other words, harmonic coordinates by themselves are 
not following the overall rotation pattern of the black- 
hole binary system. 

Binary systems provide a good example of almost- 
stationary configurations. Take for instance the well 



known pulsar 1913+16: it would be a perfect clock if 
we could just neglect the (very small) energy loss due 
to gravitational radiation, getting then a sort of steady 
system. One would like to choose a co-rotating frame in 
order to get a clean view of symmetry deviations. This 
case is representative of many other situations in which 
there is not a exact symmetry. In these cases, the idea of 
approximate symmetry, or that of almost-Killing vectors, 
would be of great help in Numerical Relativity applica- 
tions. 

A precise implementation of the concept of almost- 
symmetry has been provided by Matzner |ll| . Start- 
ing from a variational principle, it defines a measure of 
the symmetry deviation of any given spacetime. This 
idea has been applied by Isaacson to the study of 
high-frequency gravitational waves, by defining a steady 
coordinate system in which the radiation effects can be 
easily separated from the background metric. More re- 
cently [Kj , the same measure has been considered as an 
inhomogeneity index of the spacetime, which can be re- 
lated with some entropy concept. 

We are not interested here, however, in studying the 
spacetime properties or in comparing different space- 
times. We will focus instead in characterizing motions in 
arbitrary spacetimes. By a motion we mean a congruence 
of time lines, that can then be associated to the world 
lines of a system of observers. For instance, geodesic 
motions (associated with freely falling observers) or har- 
monic motions (associated with the observers at rest in 
harmonic coordinate systems). Our goal is then to char- 
acterize the motions that correspond to the physical idea 
of almost-symmetry and to study the adapted coordinate 
systems with a view to Numerical Relativity applications. 

In this sense, we will see that our approach is more 
directly related with the almost-Killing equation. This is 
a generalization of the Yano-Bochner equation 01 



v„ + ] = , 



(3) 



which has been considered by York fig and others in or- 
der to identify physically meaningful tensor components 
in asymptotically flat spacetimes (transverse-traceless 
decomposition). The same idea has also been applied to 
the identification of asymptotic Killing vectors in Kerr 
spacetime [l6| . 
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These Killing-like equations are briefly reviewed in the 
next section, where the almost-Killing equation is de- 
rived from a variational approach. The main new re- 
sults start however in the third section, where harmonic 
almost-Killing motions are shown to provide a convenient 
generalization of the standard harmonic motions. This 
generalization is implemented through a true vector con- 
dition, in contrast with the standard harmonic case, in 
which one deals instead with the set of four scalar condi- 
tions Q. This point is illustrated by considering spheri- 
cally symmetric spacetimes, where spherical coordinates 
are incompatible with the standard harmonic condition 
but perfectly allowed by the proposed generalization. 

In section four, we consider the related problem of us- 
ing harmonic almost-stationary motions as gauge condi- 
tions, with a view to Numerical Relativity applications. 
The adapted coordinate system is used in order to show 
the hyperbolicity of the full system: Einstein's field equa- 
tions plus gauge conditions. Of course, any hyperbolicity 
proof requires a specific formulation for the field equa- 
tions. We have chosen here the Z4 formalism [l^ j us t 
for simplicity, although the proposed almost-stationary 
gauge condition should also work out with other hyper- 
bolic formalisms. The specific implementation we pro- 
vide could be used as a guide for any other particular 
choices. In order to illustrate this point, we also show 
the hyperbolicity of the standard ADM free-evolution 
approach when supplemented with the proposed almost- 
stationary gauge condition. 

II. ALMOST-KILLING VECTOR FIELDS 
A. Killing-like equations 

Killing vectors can be defined as the solutions of the 
Killing equation: 

A (g» v ) = V„ + V„ = . (4) 

Their physical meaning can be better understood by con- 
sidering an adapted coordinate system. In the timelike 
case, for instance, we can choose the time lines to be the 
integral curves of £ and the time coordinate to be the spe- 
cial choice of the affine parameter on these curves such 
that 

£ = d t • (5) 
Then, the Killing equation (JIJ reads simply 

dt V = > ( 6 ) 

meaning that the metric is stationary, so that spacetime 
geometry is preserved along the integral curves of £. 

A well known generalization of the Killing equation (@J 
is given by the Affine Killing vectors ( AKV) , namely the 
solutions of 

V P [ £ 5 (g^) ] = . (7) 



The physical meaning is again more transparent if we ex- 
press it the adapted coordinate system JSJ). Then, equa- 
tion amounts to 

dt r% - o . (8) 

This means that the affine structure of the spacetime, 
given by the connection coefficients , is preserved 

along the integral curves of £. 

A interesting subset of AKV is that of the Homothetic 
Killing vectors, defined as the solutions of 

(flV) = 2 g M „ (9) 

(the factor two in the right-hand-side can be changed 
to any non-zero value by a suitable rescaling of £). 
The physical relevance of the Homothetic Killing vectors 
comes from the invariance of equation under a rescal- 
ing of the metric. This translates the idea that there is no 
preferred length (time) scale in the spacetime, allowing 
scale-invariant (self-similar) processes to develop. These 
processes have been seen to arise in connection with crit- 
ical phenomena in General Relativity 0, 0] . 

B. A variational- principle approach: almost-Killing 
vectors 

We will consider here a further generalization of 
the almost-Killing equation (AKE) given by 

v.t^^'-^v-or]^, (io) 

where the round brackets denote symmctrization. The 
solution space includes Killing vectors and AKV for any 
value of the constant A. The simplest parameter choice 
(A = 0), corresponds to the Yano-Bochner equation [l4j . 
The case A = | corresponds to the Conformal AKE, 
which includes conformal Killing vectors (A — -| for a 
three-dimensional manifold, as in ref. |15| ). 

The term 'almost-Killing' is justified by the fact that 
the AKE equation H1(J|I can be obtained from a standard 
variational principle 

SS = , S = J L ^5 d A x , (11) 

where the Lagrangian density L is given by 

i-^ ; ^ (p;CT) -^ (v-o 2 , (12) 

and the variations of the field £ are considered in a fixed 
spacetime. The covariant conservation law l|l(Jfl provides 
then a precise meaning to the heuristic concept of ap- 
proximate Killing vectors. This was not obvious a pri- 
ori, because the Lagrangian (112[1 is not a positive-definite 
quantity: the outcome of the minimization process was 
not granted to include the zeros of i|12|) . 
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The mathematical structure of the AKE is more trans- 
parent if we rewrite it in the equivalent form 

□£ M + R^u e + (1 - A) V M (V • = , (13) 

where we have just reversed the order of some covariant 
derivatives in IjlUfl ■ This 'wave equation' form of the AKE 
can be alternatively obtained from the Lagrangian 

L' = t P -, a i p -> a ~R P A P C + (1-A) (V-?) 2 , (14) 

which is of course equivalent to the original one, modulo 
a four-divergence: 

L' = 2L - V„ [ (£ ■ V) r - r (V • ] . (15) 

The principal symbol of the differential operator in ei- 
ther form of the AKE can be written in Fourier space as 

k 2 5% + (1 - A) k^K , (16) 

so that: 

• The characteristic hypersurfaces are the light cones 
(fe 2 = 0). 

• The symbol is singular for A = 2. For vac- 
uum spacetimes, this case corresponds to Maxwell's 
equations for the electromagnetic potential [lj| . A 
supplementary condition (like the 'Lorentz condi- 
tion' V • £ = 0) would be then required in order to 
get a unique solution. 

• On non-characteristic hypersurfaces, the symbol 

can be algebraically inverted for A ^ 2. 

The last point is a strong indication of the existence of 
solutions for £ in any given spacetime, for every set of 
non-characteristic initial data. Of course, this is not a 
100% rigorous proof because the straightforward passage 
to Fourier space ignores the coordinate dependence of 
the metric (the standard 'frozen coefficients' approach). 
But this gives us a sound basis for assuming in what 
follows that solutions for £ may be constructed in any 
given spacetime. 

To be more specific, the initial value problem can be 
expressed as follows 

• We can freely choose the values of £ on a given 
initial hypersurface (let us say t = 0). The space 
derivatives of £ can then be computed from these 
values. 

• The time derivative of £ can also be freely specified 
on the initial hypersurface. In this way, we have the 
full set of first covariant derivatives at t = 0. 

• In order to propagate these values along the time 
lines, we must compute the second time derivative 
of £ from the second order equation l|13|l . This can 



always be done for the space components £j pro- 
vided that g 00 7^ 0, meaning that the initial hyper- 
surface is not tangent to the local light cone. In 
the case of the £o component, we must require in 
addition that A ^ 2, as it was to be expected from 
the results of the preceding paragraph. 

In the timelike case, the integral curves of £ can be in- 
terpreted as the world-lines of a set of observers. As far as 
there is one solution for every set of (non-characteristic) 
initial data, we can interpret the set of solutions as pro- 
viding a set of motions that minimize the deviation from 
isometry along the congruence of time lines. This justifies 
the name of 'almost-stationary motions' for the timclike 
solutions of the AKE. 



III. HARMONIC ALMOST-STATIONARY 
MOTIONS 

We will consider now the particular parameter choice 
A = 1, in which the principal part of equation i|13[) is 
harmonic, that is 

ne+R» v e = o, (17) 

so that the resulting timelike solutions will be called 
'harmonic almost-Killing motions'. The conservation-law 
version (|10|l can then we written as 

V„ [-L c t (Vs«T)] = 0. (18) 

Notice that the Lagrangian L' in this case, namely 

L' = Z P ;A P '' T -R P A P ¥ , (19) 

gets an interesting 'kinetic minus gravitational' form. 

Both forms (|T7jl and ifTHl) of the Harmonic AKE equa- 
tion (HAKE) suggest a close relationship with harmonic 
coordinates. This relationship is again more transparent 
in the adapted coordinate system, where l|18(l leads to 

g p° d t r% = (20) 

(compare with Eq. [S] for AKV), whereas the harmonic 
coordinates condition reads just 

= gr° T» pa = (21) 

in adapted coordinates. One can then write 12011 as 

d t T» = v» pa d t or) , (22) 

so that it is clear that the flow associated with harmonic 
coordinates will provide a first integral for the HAKE 
in the weak field limit (where only linear terms are re- 
tained) . This fact can be relevant for the characterization 
of the gravitational waves degrees of freedom in asymp- 
totically flat spacetimes 
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Note also that the HAKE (TSJ} is a true vector equation, 
whereas the standard harmonic condition Q is rather a 
set of four scalar equations. This difference is important, 
because the congruence of time lines in a given motion is 
defined by its tangent vector field. One can, for instance, 
relabel the particular time lines by an arbitrary time- 
independent coordinate transformation, namely 

x l = f(y r ) 1,7 = 1,2,3, (23) 

while keeping the same expression JJJJ for the timelike 
tangent vector. The transformation l|23[) allows one to 
select the type of space coordinate system (cylindrical, 
spherical, or whatsoever) which is more adapted to any 
specific problem. In numerical simulations, the transfor- 
mation (|23J) corresponds to the freedom of choosing an 
arbitrary space coordinate system on the initial hyper- 
surface. 

This is not the case for the standard harmonic coordi- 
nates choice. In order to illustrate this, we will consider 
for instance a spherically symmetric line element, namely 

ds 2 = -a 2 dt 2 + X 2 dr 2 + Y 2 [d6 2 + sin 2 {9) dip 2 ] , (24) 

where all the metric functions (a, X, Y) depend only on 
(i, r). In this case, the time and radial components com- 
ponents of the HAKE (|20[l provide conditions for the cor- 
responding (i, r) coordinates, whereas the angular com- 
ponents are identically satisfied in the adapted coordi- 
nate system form <|20[) . In the standard harmonic case, 
however, one gets 

r- = o, r* = -^^o, (25) 

so that spherical coordinates happen to be incompatible 
with the harmonic condition 1|21|1 . 

This lack of versatility of the standard harmonic co- 
ordinates can be a serious drawback in Numerical Rel- 
ativity applications, where one could be unable to fully 
adapt the coordinate frame to the features of the physical 
system under consideration. We hope that the proposed 
almost-stationary generalization l|20|l will contribute to 
avoid this complication. 



IV. ALMOST-STATIONARY GAUGE 
CONDITIONS 

Standard harmonic motions, as defined by 1|21|) . have 
been used recently in advanced Numerical Relativity ap- 
plications 0, S IE 0HE3- Note however that in this 
context one is building up the coordinate system and the 
spacetime itself at the same time, whereas the spacetime 
was supposed to be given in the previous sections. To 
be more precise, Einstein's field equations must be 
coupled with the gauge condition. In the harmonic case, 
this condition is given by eq. 1|21[1 . 



The principal part of Einstein's equations can be writ- 
ten in the DeDonder-Fock form [lj, Q , namely 

□ 9^ ~ - d v T^ = ■■■ (26) 

so that, by choosing the spacetime coordinates to be the 
solutions of , one ensures the vanishing of and the 
field equations can be relaxed to a system whose principal 
part consists in a scalar wave equation for every metric 
component, namely 

□ 9ftu = ... (27) 

Note that the metric in 1(271 |27|) is over-determined, 
because one gets in all 14 equations for only 10 metric 
components. This can be better understood by introduc- 
ing a 'zero vector' as an additional dynamical field (Z4 
system ^7|)> so that the field equations read 

R„u + V„Z U + V„Z M = 8 7T (r M „ - - T 9liv ) . (28) 

The principal part of the Z4 field equations in the 
DeDonder-Fock form reads now 

□ 9^ - d^Tv + 2Z U )- d„(r M +22,) = - , (29) 

so that the relaxed system (|27fl is recovered by setting 
the values of Z 11 to be 

z^ 1 = -i r", (30) 

which can be considered just an extension of the har- 
monic coordinates condition lJT]l. namely 

□ $^=2^. (31) 

As far as one recovers in this way exactly the same 
relaxed system (I27|l . one gets exactly the same solutions 
for the metric. The extra quantities Z^ just allow us to 
monitor to which extent the harmonic coordinates con- 
dition (|30|l . when considered as a constraint on the com- 
puted metric, is actually verified. In this sense, it is useful 
to consider the four-divergence of the Z4 field equations 
(|2*5jl . namely 

□ Z ^ + Z v = , (32) 

which can be interpreted as the constraint-propagation 
law. It follows that the constraint-violation vector Z^ 
obeys precisely the HAKE l(T7|) . 

Note that the Z4 system has been used here just as 
a convenient analysis tool in order to discuss the over- 
lap between the 10 equations of the relaxed system and 
the 4 equations of the harmonic coordinates condition. 
The same overlap will occur when replacing the stan- 
dard harmonic coordinates condition H21[) by the HAKE 
(J20J- This suggests to consider, in the Z4 context, the 
analogous replacement of l|30|l by 

r3 i (rM+2 9 t Z" = 0. (33) 



5 



A covariant expression for l|33(l is given by 

V„ [± A (V5<?H]=2A • (34) 
The principal part of (|33[1 can be written now simply 

as 

d t (r" +2Z") = ... , (35) 

so that (the principal part of) the full system (J2SIE5J) gets 
a triangular form. The characteristic lines can easily be 
identified: 

• the time lines, as it follows from the fact that the 
'Gamma sector' equations ()35[) are yet in diagonal 
form 

• the light rays, as it follows from the fact that the 
only non-diagonal terms in the 'metric sector' equa- 
tions H29f) are just coupling terms with the Gamma 
sector, which is itself in diagonal form. 

We can conclude that the full differential system 
formed by the Z4 system igSJ and the (extended) HAKE 
condition l|33|) is hyperbolic. The only trouble can 
arise at the specific points where the harmonic almost- 
stationary vector £ is light-like (£ 2 = 0) , so that the cor- 
responding time line gets tangent to the light cone. The 
non-diagonal coupling terms in l|29[) would then prevent 
the full diagonalization of the characteristic matrix at 
this specific point. This would be for instance the case 
of the apparent horizon in stationary spacetimes, when 
the Killing vector is selected as a solution of the HAKE 
equation. 

Let us note again that we are using here the Z4 formal- 
ism just as an analysis tool, which allows us to monitor 
the evolution of the constraint violations. The proposed 
coordinate conditions l|2U|) are actually independent of 
the hyperbolic formalism one likes to choose for the field 



equations. In order to illustrate this point, let us take for 
instance the standard ADM free-evolution approach: 

• The original field equations (|26[1 are considered. 
However, only the space components 

□ 9ij - diTj - djTi = • • • (36) 

are kept, because the four remaining combinations 
provide just constraints which are not solved in the 
free evolution approach. 

• The almost-stationary conditions (|2(J|) , with princi- 
pal part 

d t T° = ■■■ d t Ti = • • ■ (37) 

which provide the missing evolution equations for 
the remaining metric components g 00 , goi, respec- 
tively. 

We can see by inspection that the principal part of the 
full system (|36l 1371) is also in triangular form. The char- 
acteristic lines are again either the time lines ('Gamma 
sector') and the light cones ('metric sector'). The result- 
ing ADM system, when supplemented with the almost- 
stationary gauge condition, is then hyperbolic, provided 
that the time lines do not get tangent to the light cones. 

We can conclude that the same qualitative behavior 
is obtained both in the Z4 and in the ADM frameworks 
when using the almost-stationary gauge condition. We 
do not expect then any essential difficulty in adapting 
this coordinate condition to other hyperbolic formalisms 
which are being used in Numerical Relativity applica- 
tions. 
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